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ABSTRACT 

We present a simulation of the cosmic evolution of the atomic and molecular phases of the cold 
hy drogen gas in about 3 - 10^ ga laxies, obtained by post-processing the virtual galaxy ca talog produced 
bv | De Lucia &: BlaizotI ()2007f ) on the Millennium Simulation of cosmic structure (jSpringel et al.l 
I2005D . Our method uses a set of physical prescriptions to assign neutral atomic hydrogen (HI) and 
molecular hydrogen (H2) to galaxies, based on their total cold gas masses and a few additional 
galaxy properties. These prescriptions are specially designed for large cosmological simulations, 
where, given current computational limitations, individual galaxies can only be represented by 
simplistic model-objects with a few global properties. Our recipes allow us to (i) split total cold 
gas masses between HI, H2, and Helium, (ii) assign realistic sizes to both the HI- and H2-disks, 
and (iii) evaluate the corresponding velocity profiles and shapes of the characteristic radio emis- 
sion lines. The results presented in this paper include the local HI- and H2-mass functions, the 
CO-luminosity function, the cold gas mass-diameter relation, and the TuUy-Fisher relation (TFR), 
which all match recent observational data from the local Universe. We also present high-redshift 
predictions of cold gas diameters and the TFR, both of which appear to evolve markedly with redshift. 

Subject headings: ISM: atoms - ISM: molecules - ISM: clouds - radio lines: galaxies 



1. INTRODUCTION 

Observations of gas in galaxies play a vital role in many 
fields of astrophysics and cosmology. Detailed studies of 
atomic and molecular material now possible in the local 
Universe with radio and millimeter telescopes will, over 
the coming decades, be extended to high redshifts as new 
facilities come on line. 

Firstly, hydrogen is the prime fuel for galaxies, when 
it condenses from the hot ionized halo onto the galactic 
disks. The fresh interstellar medium (ISM) thus acquired 
mainly consists of atomic hydrogen (HI), but in particu- 
larly dense regions, called molecular clouds, it can further 
combine to molecular hydrogen (H2). Only inside these 
clouds can new stars form. Mapping HI and H2 in indi- 
vidual galaxies therefore represents a key tool for under- 
standing their growth and evolution. Secondly, the char- 
acteristic Hl-radio line permits the measurement of the 
radial velocity and velocity dispersion of the ISM with 
great accuracy, thereby leading to solid conclusions about 
galaxy dynamics and matter density profiles. Thirdly, 
and particularly with regard to next-generation radio fa- 
cilities, surveys of HI are also discussed as a powerful 
tool for investigating the large scale structure of the Uni- 
verse out to high redshifts. While such large scale sur- 
veys are currently dominated by the optical and higher 
frequency bands [e.g. S pitzer (iFang et al.l l2005l) , SPS S 
(jEisenstein et al.l l2005( l. DEEP2 ()Davis et all l2003f ). 



2dFG RS (iCole et all 120051) GALEX (|Milliard et al.1 
[2007I ). Chandra (|Gini et all [20031 )] . they may well be 
overtaken by future radio arrays, such as the S quare Kilo- 
metre Array fSKA. ICarilh fc RawHngil2004D . The lat- 
ter features unprecedented sensitivity and survey speed 
characteristics regarding the Hl-line and could allow the 
construction of a three-dimensional map of ~ 10^ HI- 
galaxies in just a few years survey time. The cosmic 
structure hence revealed, specifically the baryon acous- 
tic oscillations (BAOs) manifest in the power spectrum, 
will, for example, constrain the equation of state of 
dark energy a n order of magnitude better than possi- 
ble n owadays (jAbdalla fc Rawlings 2005; Abdalla ct al] 
l2009f ). Fourthly, deep low frequency detections will pre- 
suma bly reveal HI in the intergalactic space of the dark 
ages (iCarilli et al.l [200I - one of the ultimate jigsaw 
pieces concatenating the radiation dominated early Uni- 
verse with the matter dominated star-forming Universe. 

Typically, HI- and H2-observations are considered part 
of radio and millimeter astronomy, as they rely on the 
characteristic radio line of HI at a rest-frame frequency 
oii^ = 1.42 GHz and several carbon monoxide (CO) radio 
lines, indirectly tracing H2-regions, in the ly ~ 10^ — 10'^ 
GHz band. Such line detections will soon undergo a revo- 
lution with the advent of new radio facilities such as the 
SKA and the Atacama Large Millimeter/submillimeter 
Array (ALMA). These observational advances regarding 
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HI and H2 premise equally powerful theoretical predic- 
tions for both the optimal design of the planned facilities 
and for the unbiased analysis of future detections. 

This is the first paper in a series of papers aiming at 
predicting basic HI- and H2-properties in a large sam- 
ple of evolving galaxies. Here, we introduce a suite of 
tools to assign HI- and H2-properties, such as masses, 
disk sizes, and velocity profiles, to simulated galaxies. 
These tools are subsequently applied to the ~ 3- 10^ sim- 
ula ted evo lving galaxies in th e galaxy-catalog produced 
by iDe Lucia fc Blaizoti (j2007f ) (hereafter the "DeLucia- 
catal og" ) for the Millenni um Simulation of cosmic struc- 
ture (jSpringel et al.l[2005l ). In forthcoming pubhcations, 
we will specifically investigate the cosmic evolution of HI- 
and H2-masses and -surface densities predicted by this 
simulation, and we will produce mock-observing cones, 
from which predictions for the SKA and the ALMA will 
be derived. 

Section [2] provides background information about the 
DeLucia-catalog. In particular, we highlight the hy- 
brid simulation scheme that separates structure forma- 
tion from galaxy evolution, and discuss the accuracy of 
the cold gas masses of the DeLucia-catalog. In section [3l 
we derive an analytic model for the H2/HI-ratio in galax- 
ies in order to split the hydrogen of the cold gas of the 
DeLucia-catalog between HI and H2. We compare the 
resulting mass functions (MFs) and the CO-luminosity 
function (LF) with recent observations. Sections S] and [5] 
explain our model to assign diameters and velocity pro- 
files to HI- and H2-disks. The simulation results are 
compared to observations from the local Universe and 
high-redshift predictions are presented. In Section [6l we 
discuss some consistency aspects and limitations of the 
approaches taken in this paper. Section [7] concludes the 
paper with a brief summary and outlook. 

2. BACKGROUND: SIMULATED GALAXY 
CATALOG 

iV-body simulations of cold dark matter (CDM) on 
supra-galactic scales proved to be a powerful tool to an- 
alyze the non-linea. r evolution of cosmic structure (e.g. 
ISpringel et al.ll2006[ ). Starting with small primordial per- 
turbations in an otherwise homogeneous part of a model- 
universe, such simulations can quantitatively reproduce 
the large-scale structures observed in the real Universe, 
such as galaxy clusters, filaments, and voids. These sim- 
ulations further demonstrate that most dark matter ag- 
gregations, especially the self-bound haloes, grow hier- 
archically, that is through successive mergers of smaller 
progenitors. Hence, each halo at a given cosmic time 
can be ascribed a "merger tree" containing all its pro- 
genitors. One of the most prominent sim ulations is 
the "Millennium run" (jSpringel et al.l [20051 ). which fol- 
lowed the evolution of 2160^ « 10^° particles of mass 
8.6 • 10*Mq//i over a redshift range z = 127 ^ in 
a cubic volume of (500Mpc//i)^ with periodic bound- 
ary conditions. The dimensionless Hubble parameter h, 
defined as Hq = 100 /ikms^^ Mpc^^, was set equal to 
h — 0.73, and the other cosmological parameters were 

chosen as f^mattcr = 0.25, fibaryon = 0.045, f^A = 0.75, 

0-8 = 0.9. 

Despite impressive results, modern A'^-body sim- 
ulations of CDM in comoving volumes of order 
(500Mpc/ft,)^ cannot simultaneously evolve the detailed 



substructure of individual galaxies. The reasons are 
computational limitations, which restrict both the mass- 
resolution and the degree to which baryonic and radia- 
tive physics can be implemented. Nevertheless, an ef- 
ficient approximate solution for the cosmic evolution of 
galaxies can be achieved by using a hybrid model that 
separates CDM-dominated structure growth from more 
complex baryonic physics (jKauffmann et all 1 19991 ). The 
idea is to first perform a purely gravitational large-scale 
A^-body simulation of CDM and to reduce the evolving 
data cube to a set of halo merger trees. These dark mat- 
ter merger trees are assumed independent of the bary- 
onic and radiative physics taking place on smaller scales, 
but they constitute the mass skeleton for the forma- 
tion and evolution of galaxies. As a second step, each 
merger tree is populated with a list of galaxies, which 
are represented by simplistic model-objects with a few 
global properties (stellar mass, gas mass, Hubble type, 
star formation rate, etc.). The galaxies are formed and 
evolved according to a set of physical prescriptions, often 
of a "semi-analytic" nature, meaning that galaxy proper- 
ties evolve analytically unless a merger occurs. This hy- 
brid approach tremendously reduces the computational 
requirements compared to hydro-gravitational A'^-body 
simulations of each galaxy. 

Croton et al. (2006) were the first to apply this hybrid 
scheme to the Millennium Simulation, thus producing a 
catalog with ~ 1.1 • 10^ galaxies in 64 time steps, corre- 
sponding to about ~ 3 • 10^ evolving galaxies. Th i s cata - 
log was further improved bv lDe Lucia fc Blaizoti (|2007| ). 
giving rise to the DeLucia-catalog used in this paper. 
The underlying semi-analytic prescriptions to form and 
evolve galaxies account for the most important mecha- 
nisms known today. In brief, the hot gas associated with 
the parent halo is converted into galactic cold gas accord- 
ing to a cooling rate that scales with redshift and depth 
of the halo potential. Stars form at a rate proportional to 
the excess of the cold gas density above a critical density, 
below which star formation is suppressed. In return, su- 
pernovae reheat some fraction of the cold gas, and, if the 
energy injected by supernovae is large enough, their ma- 
terial can escape from the galaxy and later be reincorpo- 
rated into the hot gas. In addition, when galaxies become 
massive enough, cooling gas can be reheated via feed- 
back from active galaxy nuclei (AGNs) associated with 
continuous or merger-based black hole mass accretion. 
These basic mechanisms are completed with additional 
prescriptions regarding merger-related starbursts, mor- 
phology changes, metal enrichment, dust evolution and 
change of photometric pro perties (see lCroton et al.ll2006t 
iDe Lucia fc Blaizoti l2007l ). The free parameters in this 
model were adjusted such that the simulated galaxies at 
redshift z — &t the joint luminosity /color / morphology 
distribution of observed low-redshift galaxies |Cole et al.l 
I2OOII: iHuang et al.ll2003L iNorberg et alll200l ). A good 
first order accuracy of the model is suggested by its abil- 
ity to rep roduce the observe d bulge-to-black hole mass 
relation dHaring &: Ri^l2004[ ). the TuUy-Fisher relation 
(jGiovanelli et al.lll997D . a nd the cold gas meta llicity as 
a function of st ellar mass ([Tremonti et al.ll2004 see also 
Figs. 4 and 6 in lCroton et al.ll2006l ). 

Some galaxies in the DeLucia-catalog have no corre- 
sponding halo in the Millennium Simulation. Such ob- 
jects can form during a halo merger, where the resulting 
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halo is entirely ascribed to the most massive progenitor 
galaxy. In the model, the other galaxies continue to ex- 
ist as "satellite galaxies" without haloes. These galaxies 
are identified as "type 2" objects in the DeLucia-catalog. 
If the halo properties of a satellite galaxy are required, 
they must be extrapolated from the original halo of the 
galaxy or estimated from the baryonic properties of the 
galaxy 

We emphasize that the semi-analytic recipes of the 
DeLucia-catalog are simplistic and may require an ex- 
tension or readjustment, when new observational data 
become available. In particular, r ecent observations 
(|Bigiel et all 120081 : iLerov et all l2008l ) suggest that star 
formation laws based on a surface density threshold are 
suspect, especially in low surface density systems. More- 
over, galaxies in the DeLucia-catalog with stellar masses 
A^stars below 4 • 10^ Mq typically sit at the centers of 
haloes with less than 100 particles, whose merging his- 
tory could only be followed over a few discrete cosmic 
time steps. It is likely that the p hysical properties o f 
these galaxies are not yet converged. ICroton et al.l ()2006f ) 
noted that especially the morphology (colors and bulge 
mass) of galaxies with Mgtars ^ 4 • 10^ Mq is poorly 
resolved, since, according to the model, the bulge for- 
mation directly relies on the galaxies' merging history 
and disk instabilities. Nevertheless, the simulated cos- 
mic space densities of stars in early-type and late-type 
galaxies at redshift z = and the cosmic star formation 
history are c onsistent with obse rvations (see figures and 
references in ICroton et al.|[2006f) . It is therefore probable 
that at least the more massive galaxies in the simulation 
are not significantly affected by the mass-resolution and 
the simplistic law for star formation. 

In this paper, we post-process the DeLucia-catalog to 
estimate realistic HI- and H2-properties for each galaxy. 
Our prescriptions will make use of the cold gas masses 
given for each galaxy in the DeLucia-catalog, and hence 
it is crucial to verify these cold gas masses against 
current observations. Based on a new estimation of 
the H2-MF, we have recently calculated the normal- 
ized density of cold neutral gas in the local Universe as 
rj°b| ^ (4.4 ± 0.8) . 10-4 fe-i, hence »g^| 6.0 • IQ-'^ 

for h = 0.73 (|Obreschkow fc Rawlingj|2009f ). This value 
was obtained by integrating the best fitting Schechter 
functions of the local HI-MF and H2-MF and therefore 
it includes an extrapolation towards masses below the 
respective detection limits of HI and H2 . The simulated 
local cold gas density fig™ of the DeLucia-catalog, ob- 
tained from the sum of the cold gas masses of all galaxies 
at redshift z = 0, exceeds the observed value by a factor 
( = = 1.45, as shown in Table[U For compar- 

ison, this table also lists the cold gas densities of other 
galaxy catalogs produced for the Mil lennium Simulation , 
using different semi-analytic recipes (jBertone et al.|[2007l) 
and different sche mes for the con struction of dark matter 
merger trees (|Bower et al.ll2006[ ). 

There are plausible reasons for the excess of cold gas 
in the DeLucia-catalog compared to observations. Most 
importantly, the semi-analytic recipes only distinguish 
between two gas phases: the hot (T « 10^ — 10 K) 
and ionized material located in the halo of the galaxy 
or group of galaxies, and the cold (T « 10^ — 10^ K) 
gas in galactic disks. However, recent observations have 



Catalog 




Qsim /Qobs 
^ 'sas / ^ 'sas 


De Lucia & Blaizot f2007~) 


8.7 ■ 10~* 
14.8 ■ 10-* 
9.0 ■ 10-* 


1.45 
2.45 
1.50 


Bower et al. (2006) 


Bertone et al. (20071 



TABLE 1 

Normalized cold gas densities at 2 = of three different 

semi-analytic galaxy simulations applied to the 
Millennium Simulation of cosmic structure. The rightmost 

COLUMN shows THE MULTIPLICATIVE OFFSET FROM THE OBSERVED 
VALUE AS DETERMINED BY ObRESCH KOW fc RAW LINGS (2009). 

clearly revealed that some hydrogen in the disk of the 
Milky W ay is warm (T 10'* K) and ionized, too. For 
example . iRevnoldi ()2004f ) analyzed faint optical emission 
lines from hydrogen, helium, and trace atoms, leading to 
the conclusion that about 1/3 of all the hydrogen gas in 
the Local Interstellar Cloud (LIC) is ionized. If this were 
true for all the gas in disk galaxies, one would expect a 
correction factor around 1.5 between simulated disk gas 
and cold neutral gas. Justified by this considerations, we 
decided to divide all the cold gas masses in the DeLucia- 
catalog M^^^'^"'^ by the constant C = 1-45 in order to 
obtain more realistic estimates, 



71 ,r — /^^l TixDcLucia 
JWgas — I, -"'^gas 

3. GAS MASSES AND MASS FUNCTIONS 



(1) 



In this section we establish a physical prescription to 
subdivide the cold and neutral hydrogen mass A/jj = 
A^Hi+A^H2 of a galaxy into its atomic (HI) and molecular 
(H2) component based on the observed and theoretically 
confirmed relation between lo cal gas pressure and local 
molecular fraction (lE hncRrc en|[l993t iBlitz fc: Rosolowskvl 
I2OO6I: ILerov et al.|[2008t iKrumholz et al.ll2009[ ). This pre- 
scription shall be applied to the DeLucia-catalog. The re- 
sulting simulated HI- and H2-mass functions (MFs) and 
the related CO-luminosity function (LF) will be com- 
pared to observations in the local Universe. 

3.1. Prescription for subdividing cold gas 

Before addressing the sub-composition of cold hydro- 
gen, we note that the total cold hydrogen mass Mh can 
be inferred from the total cold gas mass Mgas by a con- 
stant factor A/h — 0.74A'/gas, which corresponds to the 
universal abundance of hydrogen (e.g. lArnett 1996 1 that 
changes insignificantly with cosmic time. The remain- 
ing gas is composed of helium (He) and a minor frac- 
tion of heavier elements, collectively referred to as met- 
als (Z). The DeLucia-catalog gives an estimate for the 
metal mass in cold gas Mz, and hence we shall compute 
the masses of cold hydrogen and He as 



Af H = 

AfHo = 



(A'/gas 
(Afgas 



Mz) • /3, 
A/z) • (1 



(2) 



where the hydrogen fraction /3 ~ 0.75 is chosen slightly 
above 0.74 to account for the subtraction of the 1-2% 
metals in Eqs. 

The subdivision of the cold hydrogen mass A/jj = 
A/hi + Mb_2 depends on the galaxy and evolves with cos- 
mic time. We shall tackle this complexity using the vari- 
able Hz/ffl-ratio Rf^^'^''^ = MhJMhi, hence 



A/hi - A/h-(1 
A/h, = A/h-(1 



.galaxy N_i 
mol / ' 
.galaxy 
mol / 



(3) 
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Detailed observations of HI and CO in nearby regu- 
lar spiral galaxies revealed that virtually all cold gas of 
these galaxies resides in fiat, often approximately axially 
symm etric, disks (e.g. IWalter et al.l l2008t iLerov et"an 
|2008| 1. CO- maps rec ently obtained for five nearby ellipti- 
cal galaxies ()Youneli2002. ) show that even these galaxies, 
who carry most of their stars are in a spheroid, have most 
of their cold gas in a disk. There is also empirical evi- 
dence, that most cold gas in high- redshift galaxies resides 
in disks (e.g. iTacconi et al.l 120061. Based on these find- 
ings, we assume that galaxies generally carry their cold 
atomic and molecular cold gas in flat disks with axially 
symmetric surface density profiles Shi(?') and Eh2(''), 
where r denotes the galactocentric radius in the plane of 
the disk. Using these functions, Rf^li^^ can be expressed 



R'. 



galaxy 
mol 



27r/p°°drrSH,(r) 
27r drrT,m{r) 



(4) 



To solve Eq. ([4]), we shall now derive an analytic model 
for EHi(r) and Sh2('')- To this end, we analyzed the ob- 
served dens i ty pro files T,ui{r) and T,H2{r) presented by 
iLerov et al.l (|2008f l for 12 nearby spiral galaxies of The HI 
Nearby Galaxy Survey (THINGS) ^ In general, the sur- 
face density of the total hydrogen component (HI+H2) 
is well fitted by a single exponential profile. 



T,m{r) + SH2(r) = Eh exp(-r/rdisk), 



(5) 



where rdisk is a scale length and Sr = / {^nr'^^^^.) 
is a normalization factor, which can be interpreted the 
maximal surface density of the cold hydrogen disk. 

Eq. (m can be solved for T,Hi{r) and EhsI'^): if we know 
the local H2/HI-ratio in the disk, i.e. the radial function 
^moi('') = SH2('^)/SHi(r). Following the theoretical pre- 
dictio n that i?mni (r) sc a les as some power of the gas pres- 
sure (|Elmegreenlll993D . iBlitz fc Rosolowskvl (|2006[ ) pre- 
sented compelling observational evidence for this power- 
law based on 14 nearby spiral galaxies of various types. 
Perhaps the most complete empirical study of Rm n Ar) to - 
day has recently been pubhshcd by ILerov et al.l (|2008D . 
who analyzed the correlations between -Rmoi(^) and var- 
ious disk properties in 23 galaxies of the THINGS cat- 
alog. This study confirmed the power-law relation be- 
tween Rmoiir) and pr essure. On theoretical grounds, 
iKrumholz et al.l (|2009( ) argued that Rmoi{r) is most fun- 
damentally driven by density rather than pressure. How- 
ever, by virtue of the thermodynamic relation between 
pressure and density, it is, in the context of this paper, 
irrelevant which q uantity is considered, a nd the density- 
law for Rmoiir) bv IKrumholz et al] ( 2009f) is indeed con- 
sistent wit h the p ressure-laws by iBlitz fc Rosolowskvl 
(|2006D and ILerov et al. (20081 ). Here, we shall apply the 
pressure-law 



i?moi(r) = [P{r)/P, 



(6) 



where P{r) is the kinematic midplane pressure outside 
molecular clouds, and P, — 2.35 • ip-is g^j^ j ^ q, = 93 
are empirical values adopted from ILerov et al.l (|2008l ). 



^ In total ILerov et al.l II2008I ) analyzed 23 galaxies. Here, we 
only use the 12 galaxies, for which radial density profiles are 
provided for both HI and H2 (based on CO(2-l) or CO(l-O) 
measurements), an d we subtract the Helium-fraction included by 
ILerov et al] 1 I200I) . 



lElmegreenI (|1989t ) showed that the equations of hy- 
drostatic equilibrium for an infinite thin disk with gas 
and stars exhibit a simple approximate solution for the 
macroscopic kinematic midplane-pressure P{r) of the 
ISM, 



s{r){i:,Ur)+Ur)j:'^lt{r)), (7) 

where G is the gravitational constant, I]gas(?') is 
the surface density of the total cold gas component 
(HI-|-H2-|-He-|- metals), ^^^^^^(r) is the surface density 
of stars in the disk (thus excluding the bulge stars of 
early- type spiral galaxies and elliptical galaxies), and 
fa{r) = CTgas/cTstarsjZ is the ratio between the vertical 
velocity dispersions of gas and stars. The impact of su- 
pernovae and other small-scale effects on the gas pres- 
sure are implicitly included in Eq. ([7]) via the veloc- 
ity dispersion CTgas- For ^stais — 0; Eq. (O reduces to 
P{r) = 0.5 7rG'Sgas(?')^, which is sometimes used as an 
appr oximation for the I SM pressure in gas-rich galaxies 
(e.g. 'Crosthwaite fc Tur ner 2007). 
To simplify Eq. ([7]), wc note that Sgas(»') can be 

expressed as Sgas(?') = M^as/(27rr^i^^) exp(-r/rdisk), 
which is identical to Eq. ^ up to the constant fac- 
tor correcting for helium and metals. To find a simi- 
lar expression for I]gj|^g(r), we analyzed the stellar sur- 
face densities Estars(?') of the 12 THINGS galaxies men- 
tioned before. In agreemen t with many other stud- 
ies (e.g. ICourteau et al]|1996l ). we found that I]stars(?') 
is generally well approximated by a double exponen- 
tial profile, i.e. the sum of an exponential profile Sg^ara'' 
for the bulge and an exponential profile Sgtars the 
disk. On average, the scale length of the stellar disk 
fdisk is 30%-50% smaller than the gas scale length 
''disk, which traces the fact that stars form in the more 
central H2-dominated parts of galaxies. Indeed, sev- 
eral observational studies revealed that the stellar scale 
length is nearly identical to that of molecular ga s (e.g. 
lYoung et al.lll995t iRegan et al.ll2001l : ILerov et aLl HoOB). 
and hence smaller than the scale length of HI or the scale 
length rdisk of the total cold gas component. For simplic- 
ity, we shall here assume rdisk = 2fdisk for all galaxies, 
such that Sdi|k^(r) = 4 Mstars/(27r r^^J exp(-2 r/rdisk). 
Finally, we approximate the dispersion ratio /^(r) 
as /cr(r) = exp(r/rdisk), where f° is a constant. 
This approximation is motivated by empirical evidence 
that the gas dispersion Cgas rem ains approximately 
constant across galactic disc s (e.g. iDickey et aO 119901 : 
iBoulanger fc Viallefondl Il992t ILerov et al.ll2008f ). com- 
bined with theoretical and observational studies showing 
that the stellar velocity dispersion CTstars.z decreases ap- 
proximately exponentially with a scale length t wice that 
of the stellar surface density (e.g. lBottemalll993[ ). Within 
those approximations Eq. ([7]) reduces to 



P(r) 



gas 



S^^disk 



^(/.)M,tt exp(-2r/rdisk), (8) 



where (fa-) = 4/^ is a constant, which can be interpreted 
as the average value of (r) weighted by the stellar sur- 

face density, since / 2^ r E,trs(0/.(0A4t's = 4/°. 
Substituting F(r) in Eq. §1 for Eq. we obtain 



SH2W ^ 
SHi(r) 



Rmoi{r) = i?moi exp(-1.6r/rdisk) (9) 
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with 



DC _ 



disk 



0.8 



(10) 



where K = G/(87rP*) = 11.3 m^kg-^. Eq. ^ reveals 
that the H2/HI-ratio i?moi(r) is described by an exponen- 
tial profile with scale length rdisk/1-6. It should be em- 
phasized that the central value R'^^^i does not necessarily 
correspond to the H2/HI-ratio measured at the center of 
real galaxies due to an additional H2-enrichment caused 
by the central stellar bulge. However, R'^^i represents 
the extrapolated cental H2/HI-ratio of the exponential 
profile, which approximates Rmoiir) in the outer, disk- 
dominated galaxy parts. 

We can now solve Eqs. (0 for the atomic and molec- 
ular surface density profiles, i.e. 



SHi(r-) : 



SH.(r): 



Sh exp(-r/rdisk) 
1 + -Rmoi'=xp("l-6^Adisk) ' 

^H-Rmoi exp(-2.6r/rdisk) 



(11) 



(12) 



1 + -^moi exp(-1.6r/rdisk) 
These model-profiles can be checked against the ob- 
served HI- and H2-density p r ofiles of the nearby galax- 
ies analyzed by ILerov et aU (|2008f) . In particular, we 
can test the limitations implied by our assumption that 
rdisk = 2fdisk- To this end, we selected two observed 
regular spiral galaxies with rdisk ~ 2fdisk (NGC 3184) 
and rdisk ~ r=disk (NGC 5505). To evaluate Eqs. pTlfT^ 
for those galaxies, we require the quantities Mg'l'f^g, Mgas, 
rdisk, and {f„). M^^f^^, M^^, and rdisk were determined 
by fitting a single exponential profile to the total cold 
gas component and a double exponential profile (bulge 
and disk) to the stellar component, a nd (/j) was chosen 
as {fa) = 0.4, i.e. the value given by lElmegreenl ()1993D 
for nearby galaxies. As shown in Fig. [1] the resulting 
model-profiles Shi {f) and {f) approximately match 
the empirical data. The fact that the fit is rather good 
for both galaxies demonstrates that the quality of the 
model-predictions does not sensibly depend on the good- 
ness of the model-assumption rdisk ~ 2fdisk- Similarly 
good fits are indee d found for most of the 12 THINGS- 
galaxies, for which ILerov et al.l (|2008D published radial 
HI- and H2-density profiles. 

Eqs. ([Tl] [T2|) can be solved for the maximal surface 
densities of HI and H2. EHi(r) exhibits its maximum at 
the radius r^f^ = 0.625 rdisk ln(3/5 • R^-^ai), as long as 
-^moi > 5/3. Galaxies in this category show an Hl-drop 
towards their center, suc h as observed in mo st galaxies 
in the THINGS catalog (jWalter et al.ll2008f ). By con- 
trast, disk galaxies with R'^^i < 5/3 have Hl-density 
profiles peaking at the center, rjjf^ = 0. Galaxies with 
such small values of R'l 



„max 
'^HI 

"•moi have low gas densities by 
virtue of Eq. (fTO]) . such as the irreg ular galaxies NGC 
4214 and NGC 3077 (see profiles in ILerov et aLll2008h . 
{r) given in Eq. (fT^ always peaks at the disk center. 
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Fig. 1. — Column density profiles of two nearby spiral galax- 
ies. Filled triangles, circles, and squares respectively r epresent the 
meas ured column density profiles of stars, HI, and H2 IILerov et al.1 
120081) . Solid lines show the best-fitting double-exponential func- 
tions for the stellar densities. Dashed lines and dash-dotted lines 
represent the predictions of our pressure-based model given in 

Eqs. imira . 

The density profiles of Eqs. ((Tn[T2)) can be substituted 
into Eq. (jl]). The exact solution of Eq. ([4]) is quite un- 
handy, but Rf^li^^ only depends on R'^^^i and an excellent 
approximation, accurate to better than 5% over the nine 
orders of magnitude = 10~^ — 10^ (covering the 

most extreme values at all redshifts), is given by 



R 



galaxy 
mol 



(3.44 i?;^ 



-0.506 



4.82 i?^ 



-1.054\-1 



(15) 



Eqs. pm [T5|) constitute a physical prescription to 
estimate the H2/HI-ratio of any regular galaxy based 
on four global quantities: the disk stellar mass Mg'j.'f^g, 
the cold gas mass Mgas, the scale radius of the cold 
gas d isk rnisk, and the d ispers ion parameter (fa-)- In 
Obre schkow fc Rawlingsl ()2009r ). we showed that the 
H2/HI-ratios inferred from this model are consistent with 
observations of nearby galaxies. Moreover, this model 
presumably extends to high redshifts, since it essentially 
relies on the fundamental relation between pressure and 
molecular fraction and on a few other physical assump- 
tions with weak or absent dependence on cosmic epoch. 
However, a critical discussion of the limitations of this 
model is presented in Sections 16.21 and 16.31 

3.2. Application to the DeLucia- catalog 

We applied the model given in Eqs. pHl [T5)) together 
with Eqs. (in [3]) to the DeLucia-catalog in order to as- 
sign HI-, H2-, and He-masses to the simulated galaxies. 
The quantities and Mz used in these Eqs. are di- 

rectly contained in the DeLucia-catalog, and Mgas was 
inferred from the given cold gas masses via the correc- 
tion of Eq. 11]). The dispersion parameter (fa-) is approx- 
imated by the constant = 0.4, consistent wit h the 
local observational data used by lElmegreenl (119891 ) (but 
see discussion in Section lO]) . 

The remaining and most subtle ingredient for our pre- 
scription of Eqs. pm [T5|) is the scale radius rdisk- This 
radius can be estimated from the virial radius rvir of the 
parent halo, but their relation is intricate. Even modern 
TV-body plus SPH simulations of galax y formation can- 
not reproduce observed disk diameters |Kaufmann et all 
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l2007f) . and hence the more simplistic semi-analytic ap- 
proaches arejhkely to require some empirical adjustment. 
iMo et ail (|1998( ) studied the case of a flat exponential 
disk in an isothermal singular halo. When assuming that 
the disk's mass can be neglected for its rotation curve, 
they find 

rdisk^ —^r^ii, (16) 

where A is the spin parameter of the halo and ^ is the 
ratio between the specific angular momentum of the disk 
(angular momentum per unit mass) and the specific an- 
gular momentum of the halo. The model behind Eq. (|16p 
does not distinguish between different scale radii for stars 
and cold gas, but it assumes a single exponential disk in 
hydrostatic equilibrium without including the effects of 
star formation. It is therefore natural to identify rdisk in 
Eq. with the cold gas scale radius rdisk of Eqs. ([TU]- 

ini). 

In the Millennium Simulation, Tvir was calculated from 
the virial mass Mvir using the relation 

Air 

M^i, ^ —rl,^. ■ 200 p,{z) (17) 

where pdz) is the critical density for closure pdz) = 
37J^(z)/(87rG) and Mvir the virial mass of the halo. For 
central haloes, Mvir was approximated as M200, i-e. the 
mass in the region with an average density equal to 
200 Pc{z)', and for sub-haloes, Mvir was approximated as 
the total mass of the gravitationally bound simulation- 
particles. 

The spin parameter A was calculated directly from the 
A^-body Millennium Simulation according to the defi- 
nition A EE Jhalo-Bhalo"^^^G"^Mhalo~^^^, whcrC Jhalo dc- 

notes the angular momentum of the halo, £'haio its en- 
ergy, and Mhaio its total mass. For the satellite galaxies 
in the DeLucia-catalog, i.e. the ones without halo (see 
Section [5]), the value of A • rvir was approximated as the 
respective value of the original galaxy halo before its dis- 
appearance. 

The only missing parameter for the calculation of rdisk 
via Eq. (|16p is the angular momentum ratio £. It is of or- 
der unity for evolved disk -galaxies (e.g. lFall fc Efstathioul 
119801: IZavala et al.l[2008h . but its exact value is uncertain 
because of the difficulty of measuring the spin of dark 
matter haloes and conv ergence issues in numer ical sim- 
ulations. For example, iKaufmann et al.l ()2007r ) showed 
that iV-body plus SPH simulations with as many as 10^ 
particles per galaxy do not reach convergence in angu- 
lar momentum, because of the difficulties to model the 
transport of angular momentum. 

Here, we shall chose ^, such that our simulation repro- 
duces the empirical relation between the galaxy baryon 
mass Mbary and the stellar scale radius rdisk ~ ?'disk/2, 
measured in the local Universe. To ensure consis- 
tency with Section 13. li we use again the data fro m the 
THINGS-galaxies analyzed by (iLerov et al.l [Mil . Of 
the 23 galaxies in this sample, we reject the 6 irreg- 
ular objects, since our model for HI and H2 assumes 
regular galaxies. The remaining 17 galaxies cover all 
spiral types and include Hl-rich and H2-rich galaxies. 
For each galaxy, we adopted the stellar scale radii fdisk 
and baryon masses Mbary = Mgtars + Mhi + A'/h2 di- 



rectly from the data presented in (jLerov et al.l l2008l )^. 
The Hubble-types T, i.e. the numerical stage indexes 
along the revised Hubble s equence of the RC2 system 
()de Vaucouleurs et al.lll976D. were dra wn from the Hy- 
perLeda database (jPaturel 'etal][20n3h . 



late-type spiral galaxies 




-1.5 -1.0 -0.5 0.0 0.5 1.0 1.5 

log(M,„,/10"'Mo) 

Fig. 2. — Relation between the baryon mass (stars+cold gas) 
and the stellar scale radius of disk galaxies. The filled and empty 
squares, respectively, represent observed early- and late-type spiral 
galaxies jLerov et al.| [2008). Typical l-tr error bars are shown in 
the bottom right corner. The empirical fit given in Eq. (1191 is 
represented by a solid line for early-type spiral galaxies (T = 2) 
and by a dashed line for late- type spiral galaxies (T = 10). The 
dotted and dash-dotted lines represent the respective relations of 
the simulated galaxies in the DeLucia-catalog, if ^ = 1. In order 
for the simulation to maxima lly align with the observations, g must 
be chosen according to Eq. I|20| l. For the latter case, 10^ random 
early- and late-type spiral galaxies of the DeLucia-catalog at 2: = 
are represented by filled and empty dots, respectively. 

The resulting empirical relation between Mbary and 
fdisk is displayed in Fig. [2l The scatter of these data 
probably underestimates the true scatter caused by all 
galaxies, since we exclusively considered non-interacting 
regular spiral galaxies. The zero-point of the mean re- 
lation between Mbary and fdisk depends on the morpho- 
logical galaxy type, as is revealed by the distinction of 
early-type and late-type spiral galaxies in Fig. [2l Given 
identical baryon masses, the stellar scale radii of early- 
type galaxies tend to be smaller than the radii of late- 
type galaxies. This trend w as also detected in other 
d ata samples (e.g. data from pCregel et all |2002| shown 
in lObreschkow fc Rawhngj|2009[ ). Several reasons could 
explain this finding: (i) early-type galaxies have more 
massive stellar bulges, which present an additional cen- 
tral potential that contracts the disk; (ii) bulges often 
form from disk instabilities, occurring preferably in sys- 
tems with relatively low angular momentum, and hence 
early-type galaxies are biased towards smaller angular 
momenta and smaller scale radii; (iii) larger bulges like 
those of lenticular and elliptical galaxies, often arise from 
galaxy mergers, which tend to reduce the specific angular 
momenta and scale radii. 

To parameterize the dependence of the scale radius 
on the morphological galaxy type, the latter shall be 

^ Stellar masses Mstars rely on 3.6/xm-maps (SINGS, 
[Kennicutt ct al.' '2003^; HI -masses use the 21 cm-maps from 
THINGS (Walter ct al. 2008), H2-m asses r ely on CO-maps (CO(2- 
1) from HERACLES. iLerov et all l2008l : CO(l-O) from BIMA 
SONG, IHelfer et al.ll200ir 
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quantified using the stellar mass fraction of the bulge, 
B = M^^^^ /Mstars- In the observed sample, the values 
of B can be approximately inferred from the Hubble type 
T. Here we shall use the relation 

6 = (10-T)V256, (18) 
which approximately parameterizes the mean behavior of 
146 moderately incli ned barred and unbarred local spiral 
galaxies analyzed bv lWeinzirl et al.l |2009). Eq. sat- 
isfies the boundary conditions B = 1 for T = —6 (i.e. pure 
spheriods) and ,B = for T = 10 (pure disks). 

We shall approximate the relation between Mbary and 
fdisk as a power-law with an additional term for the ob- 
served secondary dependence on morphological type, 

where oq, ai, and 02 are free parameters. The best fit 
to the empirical data in terms of a maximum-likelihood 
approach is given by the choice aq = 0.3, ai = 0.4, 02 = 
-0.6. 

The fit of Eq. HH) is displayed in Fig. [2] for early- 
type spiral galaxies [T = 2 ^ B = 0.25, solid fine) and 
for late-type spiral galaxies (T = 10 ^ S = 0, dashed 
line) . For comparison, the mean power-law relations for 
the early- and late-type spiral galaxies in the DeLucia- 
catalog at z = are displayed as dotted and dash-dotted 
lines for the choice ^ = 1. The simulated scale radii using 
^ = 1 are consistent with the observed ones for very 
massive galaxies (Mbary ~ 10^^ Mq), but less massive 
galaxies in the simulation turn out slightly too large if 
^ = 1. Furthermore, the morphological dependence of 
the simulation is too small compared to the observations. 
This can be corrected ad hoc by introducing a variable ^ 
that depends on both Afbary and B, i.e. 

log(0 -60 + 61 log (^^Sg^) + b2 B, (20) 

where 60, 61, and 62 are free parameters. The param- 
eters minimizing the rms-deviation between the simu- 
lated galaxies and the empirical model of Eq. ([19]) are 
60 — —0.1, 61 = 0.3, 62 = —0.6. In addition, we chose a 
lower limit for ^ equal to 0.5, in order to prevent unreal- 
istically small scale radii. 

We emphasize that Eq. (PO)) is merely an empirical cor- 
rection; this choice of ^ should not be considered as an 
estimate of the true ratio between the specific angular 
momenta of the disk and the halo, but it also ac counts for 
the imperfection of the simplistic halo model bv lMo et al.l 
(ji998), for missing physics in the semi-analytic modeling, 
and for possible systematic errors in the spin parame- 
ters A of the Millennium Simulation. The average value 
of Eq. (pn|) over all galaxies in the DeLucia-catalog is 
(^) — 0.7 (with a — 0.2), which is approximately consis- 
tent with ^ « 1 of mod ern high-resolution simulations of 
galaxy formation (e.g. iZavala et alll2008[ ). even though 
the latter still suffer from issues with the transport of 
angular momentum as mentioned above. 

Using Eqs. we estimated a scale radius rdisk 

for each galaxy in the DeLucia-catalog. A sample of 10^ 
simulated early- and late-type spiral galaxies at z = is 
shown in Fig. [51 Given rdisk as well as Mgtars, -^^gas, and 
(/^) = 0.4, we then applied Eqs. ([1 [1 [TOl [I5|) in order 
to subdivide the non- metallic cold gas mass (Mgas — Mz) 
of each galaxy into HI, H2, and He. 



3.3. Atomic and molecular mass functions 

We shall now compare the HI- and H2-masses pre- 
dicted by our model of Sections 13.11 and 13.21 to recent 
observations in the local Universe. From the view- 
point of the simulation, a fundamental output are the 
mass functions (MFs) of HI and H2, while the available 
observational counterparts a re the luminosity f unctions 
(LFs) of the Hl-emissio n line (IZwa an et al.l l2005[ ) and the 
CO(l-0)-emission line ([Keres et al...2003, ). Therefore, ei- 
ther the simulated data or the observed data need a 
luminosity-to-mass (or vice versa) conversion to compare 
the two. Section 13.31 focuses on the MFs, adopting the 
standard luminosit y-to-mass conversion for HI used by 
IZwaan et al.l (I2005D and the CO-luminosity -to-H2-mass 
conversion of lObreschkow fc Rawlingj (|2009f ) . As a com- 
plementary approach. Section 15^ will focus on the LFs, 
which will require a model for the conversion of simulated 
H2-masses into CO-luminosities. 

We define the MFs as (/)x(Afx) = dpx/dlogMx, where 
Px(Mx;) is the space density (number per comoving vol- 
ume) of galaxies containing a mass Mx of the constituent 
x (HI, H2, He, etc.). Given a mass, such as Mhi, for 
each galaxy in the DeLucia-catalog, the derivation of 
the corresponding MF only requires the counting of the 
number of sources per mass interval. We chose 60 mass 
intervals, logarithmically spaced between 10* Mq and 
10^^ M0, giving about ~ 10^ galaxies per mass interval 
in the central mass range, while keeping the mass error 
relatively small (Alog(Af) < 0.05). Since MFs combine 
units of mass and length, they generally depend on the 
Hubble constant Hq-, or the dimensionless Hubble param- 
eter h, defined as Hq = 100 /ikms~^ Mpc~^. Although 
MFs are often plotted in units making no assumption 
on h (e.g. Mq h^^ for the mass scale), this is impossible 
when observations are compared to cosmological simula- 
tions. The reason is that simulated masses in the Mil- 
lennium Simulation scale to first order as h^^, whereas 
empirical masses, when determined from electromagnetic 
fluxes, are proportional to the square of the distance and 
hence scale as h~^. For all plots in this paper we shall 
therefore use h = 0.73, which corresponds to the value 
adopted by the Millennium Simulation (see Section [2]). 

Fig. [3| displays the HI- and H2-MF of our simulation 
(solid lines), as well as the corresponding empirical MFs 
for the local Universe (points with error bars) . The 
empirical HI-MF was obtained by iZwaan et al.l (|2005l ) 
based on 4315 galaxies of the Hl-Parkes All Sky Sur- 
vey (HIPASS) a nd the empirical H2-MF was derived in 
lObreschkow fc Rawlings (2009) from the CO-lu minosity 
function (LF) presented by IKeres et al.l (|2003l ). Both 
empirical MFs approximately match the simulated data. 
We note, however, that the consistency between obser- 
vation and simulation decreases if we skip the overall 
correction of the cold gas masses in the DeLucia-catalog 
by the constant factor C, (dotted lines), which, as argued 
in Section [2| can be justified by a fraction of the disk gas 
being electronically excited or ionized. 

Our simulation slightly over-predicts the observed 
number of the largest HI- and H2-masses, i.e. the ones 
in the exponential tail of the MFs in Fig. [3l These tails 
contain the most massive systems, whose emergent lu- 
minosities are most likely to be biased by opacity and 
thermal effects. Including these effects would proba- 
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bly correct the space density of massive systems towards 
the simulated MFs. Additionally, we note that the pre- 
sented empirical MFs neglect mass measurement errors, 
which might have an important effect on the slope of 
the exponential tails. Another difference between obser- 
vations and simulation are the spurious bumps in the 
low mass range of the simulated MFs, i.e. log(MHi) ~ 
8.5 and log(MH2) ~ 8.0, where the number density 
is about doubled compared to observations. This fea- 
tur e can also be seen in the optical bj-band LF shown 
by ICroton et al.l ((2006,) and stems from an imprecision 
in the number density of the smallest galaxies in the 
DeLucia-catalog, where the mass-resolution of the Mil- 
lennium Simulation implies a poorly resolved merger his- 
tory. The over-density of sources around this resolution 
limit roughly balances the mass of even smaller galaxies, 
i.e. log(A/Hi/Mo) <C 8.0, that are missing in the simula- 
tion. 

The universal gas densities of the simulation, expressed 
relative to the critical density for closure, are = 3.4- 
10~^ and 511}™ ~ ^-^ ' lO^^i in good agreement with the 
observations fl'^^' = (3.6±0.4) ■ IQ-^ (jZwaan eFaI|[2005l) 
and fJg^; = (0.95 ±0.37) • lO'^ (lObreschkow fc Rawlingd 
[20091) . 
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Fig. 3. — Simulated galaxy MFs for HI and H2 with (solid 
lines) and without (dotted lines) the constant correction for all 
cold gas masses given in Eq. |[TJ . Filled and open squa res with error 
bars r epres ent the corresponding empirical MFs from lZwaan et ahl 
pOOSl ) and lObreschkow fc Rawlingsl l(2009h . 

Fig.dlshows our simulated HI-MF and H2-MF together 
with the MF for the cold gas metals given in the original 
DeLucia-catalog and the MF for He as trivially derived 
using Eq. ([2]). This picture reveals that in the cold gas 
of the local Universe He is probably more abundant than 
H2, but less abundant than HI. 

We shall now consider the Hl-masses in elliptical and 
spiral galaxies separately. This division is based on the 
Hubble type T, where we consider galaxies with T < 
as "ellipticals" and galaxies with T > as "spirals" - 
a rough separation that neglects other types like irregu- 
lar galaxies as well as various subclassifications. In the 
simulation, T is computed from the bulge mass fraction 
according to Eq. (fTSl) . 

The simulated HI-MFs of both elliptical and spiral 
galaxies are shown in Fig. [5l In order to determine the 
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Fig. 4. — Simulated galaxy MFs for HI (solid line), H2 (dashed 
line), cold He (dotted line), and cold gas metals (dash-dotted line). 
The HI-MF and H2-MF are identical to the solid lines in Fig. |3] 



observational counterparts, we split the HIPASS galaxy 
sample into elliptical and spiral galaxies according to 
the Hubb le types provi ded in the HyperLeda reference 
database dPaturel et al.| [2003) . For both subsa mples, the 
HI-M F was evaluated using the 1/T4iax method (jSchmidtl 
Il968l ). where T4iax was estimated from the analytic com- 
pleteness function for HIPASS, which characterizes the 
completeness of each source given its Hl-peak flux density 
Sp and integrated Hl-line flux S'int (^Zwaan ct al. 2004) . 
In order to estimate the uncertainties of the MFs, we 
derived them for 10^ random half-sized subsets of the 
HIPASS sample - a bootstrapping approach. The stan- 
dard deviation of the lO** values of \og{(f>m) for each 
mass bin was divided by V2 to estimate the l-cr errors 
of log(0Hi) for the full sample. 

Fig. [5] demonstrates that our simulation successfully 
reproduces the Hl-masses of both spiral and elliptical 
galaxies for Hl-masses greater than ~ 10^ Mq, although 
the nearly perfect match between simulation and obser- 
vation may be somewhat coincidental due to the uncer- 
tainties of the Hubble types T calculated via Eq. (fT8)) . 
For Hl-masses smaller than 10^ Mq, the morphological 
separation seems to breakdown (shaded zone in Fig. [5|). 
Indeed the Hl-mass range Afni ^ 10^ Mq approximately 
corresponds to the stellar mass range Mgtars ^ 4-10^ M©, 
for which morphology properties are poorly resolved (see 
Section 12). 

3.4. Observable HI- and CO -luminosities 

The characteristic radio line of HI stems from the hy- 
perfine energy level splitting of the hydrogen atom and 
lies at 1.42 GHz rest- frame frequency. The velocity- 
integrated luminosity of this line Lhi can be calculated 
from the Hl-mass via 



M, 



HI 



1.88 • 10^ 



^HI 



Jy km s ^ Mpc^ 



(21) 



Eq. (f2T|) neglects Hl-self absorption effects, but this is 
likely to be a proble m only for the l argest disk galax- 
ies observed edge-on ()Rao et al.lll995[ ). The strict pro- 
portionality between Lhi and Afni assumed in Eq. (|2ip 
means that the HI-LF is geometrically identical to the 
HI-MF. 
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Fig. 5. — HI-MFs for elliptical and spiral galaxies. The solid 
line and dashed line respectively represent the simulated result, 
where the galaxies have been divided in ellipticals an d spirals ac- 
cording to their Hubble type estimated using Eq. 11181 . Filled and 
open dots with error bars represent the corresponding empirical 
HI-MFs, which we derived from the HIPASS sample. The shaded 
zone represents the Hl-mass range Mjji ^ 10^ Mq, approximately 
corresponding to the simulated galaxies with poorly resolved mor- 
phologies. 

By contrast, the H2-niasses used for the empirical H2- 
MFs in Figs. [3] and [4] rely on measurements of the C0(1- 
0)-line, i.e. the 115 GHz radio line stemming from the 
fundamental rotational relaxation of the most abundant 
CO-isotopomer ^^C^^O. Here we only consider this line, 
but luminosities of other CO-lines can be estimated usin ; 
approximate emp irical line ratios (e.g. lBraine et aOlQQ 
iRighi et al.ll200l . 

The CO(l-0)-to-H2 conversion generally depends on 
the galaxy and the cosmic epoch, and it is often repre- 
sented by the dimensionless factor 



X 



_ 7Vh2 / cm 



• 10" 



(22) 



/co/(Kkms-i) 

where iVna is the column density of H2-molecules and 
/co is the integrated CO(l-0)-hne intensity per unit sur- 
face area defined via the surface brightness temperature 
in the Rayleigh- Jeans approximation. The definition of 
Eq. (1221) implies the mass luminosity relation (e.g. review 
bv .Young k. ScoviUe 19911 

A^Hs Leo 



313 



(23) 



Mq Jykms-iMpc2' 

where Leo is the velocity-integrated luminosity of the 
CO(l-O) line. 

As discussed in lObreschkow fc Rawlingsl ()2009D . the 

theoretical and observational determination of the X- 
factor is a subtle task with a long history. Most present- 
day studies assume a constant X-factor Xc, such as 



(24) 



which is typical fo r spiral galaxie s in the local Universe 
(iLerov et al.ll2008l). Bv contrast, lArimoto et all (|1996D 
an d iBoselli et al.l(|2002( i suggested that X is variable, Xv, 
and approximately inversely proportional to the metal- 
licity O/ff, i.e. the ratio between the number of oxygen 
ions and hydro gen ions in the hot ISM. U sing their data, 
we found that (jObreschkow fc Rawhng£| [2009) 

log(Xv) = (-2.9 ± 0.2) - (1.02 ± 0.05) log(0/i7). (25) 



At first sight, the empirical negative dependence of X 
on the metallicity seems to contradict the fact that 
12q16q jg optically thick for 115 GHz radiation. Indeed, 
the radiated luminosity should not depend on the den- 
sity of metals, as long as the latter is high enough for 
the r adiation to remain optically thick (jKutner fc Leun j 
19851). However, detail ed theoretical investigations (e.g. 



Malonev fc Bla^[l98l of the sizes and temperatures of 



molecular clumps were indeed able to explain, and in fact 
predict, the negative dependence of X on metallicity. 

Eq. ((251) hnks the metallicity of the hot ISM to the 
X-factor of cold molecular clouds, and it is likely a con- 
sequence of a more fundamental relation between cold 
gas metallicity and X. To uncover such a relation, we 
assume that the OjH metallicity of the cold ISM in lo- 
cal galaxies is approximated by OjH of the hot ISM. 
Given an atomic mass of 16 for Oxygen, the fact that 
hydrogen makes up a fraction 0.74 of the total baryon 
mass, and assuming that Oxygen accounts for a fraction 
of 0.4 of the ma ss of all metals (based on lArnettI 119961 : 
iKobulnickv fc Zarits kv 1999), Eq. ^ translates to 

(26) 



Xv«0.04Mgas/Mz, 



where Mgas is the total cold gas mass and Mi is the mass 



of metals in cold gas. Eq. (|26p only relates cold gas prop- 
erties to each other and therefore is more fundamental 
than Eq. (|^. 

To evaluate Eq. (pS)) for each galaxy in the simu- 
lation, we used the cold gas metal masses Mz given 
in the DeLucia-catalog. Thos e masses are reason - 
ably accurate as d e monst rated bv lDe Lucia et al.l ()2004D 
and ICroton et al.l ()2006[ ) through a comparison of the 
simulated stellar mass-metallicity relation to the em- 
pirical mass-metallicity relation obtained from 53,000 
star forming galaxies in the Sloan Digital Sky Survey 
(|Tremonti et al.|[2004f) . For most galaxies at z = the 
simulation yields metal fractions Mz/Mgas ~ 0.01 — 0.04 
in the local Universe, thus imply ing Xy » 1 — 4 in agree- 
ment with observed values (e.g. Bos elli et al.ll200^ . 

Fig. [5] displays the simulated CO-LF for the vari- 
able X-factor X^ (solid line) and the constant X-factor 
Xc (dashed hue) together with the empirical CO-LF 
(|Keres et al.ll2003[ ). adjusted to h = 0.73. The compar- 
ison supports the variable X-factor of Eq. against 
Xc — 2 (and the same conclusion is found for other con- 
stant values of Xc). Using Eq. (p6)) also has the advan- 
tage that the cosmic evolution of the X-factor due to 
the evolution of metallicity is implicitly accounted for. 
Nevertheless Eq. (|26l) may not be appropriate at high 
redshift as discussed in Section [ 



4. COLD GAS DISK SIZES 

Using the axially symmetric surface density profiles for 
HI and H2 given in Eqs. pTl [T^ . we can define the HI- 
radius rni and H2-radius of an axially symmetric 
galaxy as the radii corresponding to a detection limit 
Eq, i.e. 



SHi(mi) = So , 
Sh2(''H2) = So ■ 



(27) 
(28) 



In this paper, we chose Sq = IMqpc ^, corr esponding 
to th e deep survey of the Ursa Major group bv lVerheiienI 
(e.g. [2OOI . but any other value could be adopted. In 



10 



Obreschkow et al. 



O. T, 



o 

o 



-4 



-5 



4 5 6 7 8 

log(ico/[Jykms-' Mpc']) 

Fig. 6. — Luminosity function of CO(l-0)-emission (CO-LF) in 
the local Universe. The solid line represents the simulated CO- 
LF, obtained using the variable conversion factor Xv of Eq. I I26I I. 
and the dashed line represents the CO -LF, obtained using the con- 
stant conversion factor Xc of Eq. II24I I. Square dots and error bar s 
represent the empirical CO-LF determined bv lKeres et al.l (|2003I V 

general Eqs. ((27l [28]) do not have explicit closed-form 
solutions and must be solved numerically for each galaxy. 

Results for thi and r^^ at three epochs are displayed in 
Fig. [71 Each graph shows 10'^ simulated galaxies, drawn 
randomly from the catalog with a probability propor- 
tional to their cold gas mass. This selection rule ensures 
that rare objects at the high end of the MF are included. 
The arithmetic average of the points in each graph can 
be interpreted as the cold gas mass-weighted average of 
the displayed quantities. This average is marked in each 
graph to emphasize changes with redshift. The data in 
Fig. [7^ are shown again in Fig. [5] together with mea- 
surements of 39 spiral galaxies in the Ursa Major group 
((ycrhciien 2001). 

Figs. [7] and [5] reveal several features, which we shall 
discuss hereafter: (i) the mass-radius relation for HI is 
a nearly perfect power-law with surprisingly small scat- 
ter; (ii) in general, radii become smaller with increasing 
redshift; (iii) the evolution of the mass-radius relation is 
completely different for HI and H2. 

The first result, i.e. the strict power-law relation be- 
tween Mhi and Thi, is strongly su pported by mea sure- 
ments in the Ursa Major group ()Verheiienl 120011 see 
Fig. [5]). The best power-law fit to the simulation is 



HI 



M 



12- 







pc 



(29) 



The rms-scatter of the simulated data around Eq. 
is a = 0.03 in log-space, while the rms-scatter of the 
observations bvlVerheiienI ^200l|) is cr = 0.06. This small 
scatter is particularly surprising as the more fundamental 
relation between Mbary and fdisk shown in Fig. [2] exhibits 
a much larger scatter of a = 0.26. The square-law form 
of the power-law in Eq. ([29]) implies that the average HI- 
surface density inside the radius rui is nearly identical 
for all galaxies, which have most of their Hl-mass inside 
the radius rni. 
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Fig. 7. — Simulated mass-radius relations for HI and H2 at red- 
shifts z = 0, 4.89, 10.07, corresponding to the simulation snapshots 
63, 21, 12. Black dots represent 10^ simulated galaxies and the 
solid lines show the power-law regression for the data in Fig. [tK 
(i.e. HI at 2 = 0). The red crosses represent the cold gas mass- 
weighted averages of (MHi,r'Hl) and (Mh2 1^112) the simulation 
at each of the three redshifts. 
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Fig. 8. — Relation between Hl-mass Mjji and Hl-radius rjji for 
galaxies at redshift 2 = 0. The black dots represent 10^ simulated 
galaxies and the solid line their linear regression. The slope of this 
power-law is 0.5, thus indicating a universal average Hl-surface 
density for all disk galaxies. The simulated data are identical to 
those plotted in Fig. [Jt. Red squares show measurements in the 



Ursa Major group bv lVerheiienI 11200 Ih . who used the same defini- 
tion of rjji as this paper. 



The existence of such a constant average density of 
HI can to first order be interpreted as a consequence 
of the fact that HI transforms into H2 and stars as soon 
as its density and pressure are raised. In fact, observa- 
tions show that Ehi saturates at about 6 — 1 pc~^ 
([Blitz fc Rosolo wskv 2006|; iLerov et alJ [2008') and that 
higher cold gas densities are generally dominated by • 
Therefore, HI maintains a constant surface density dur- 



Cosmological Simulation of Neutral Hydrogen 



11 



ing the evolution of any isolated galaxy as long as enough 
HI is supplied from an external source, e.g. by cooling 
from a hot medium as assumed in the recipes of the 
DeLucia-catalog. This also explains why the power-law 
relation between Mhi and rpji remains nearly constant 
towards higher redshift in the simulation (Fig. [7K"C). 
About 1% of the simulated galaxies at redshift z — 
lie far off the power-law relation (i.e. are outside 5-cr of 
the best fit) , typically towards smaller radii (see Fig. [8|) . 
One might first expect that these objects have a higher 
Hl-surface density, while, in fact, the contrary applies. 
These galaxies have very flat Hl-profiles with most of the 
Hl-mass lying outside the radius rni, and therefore they 
would require a lower sensitivity limit than 1 Mq pc~^ for 
a useful definition of thi. Such galaxies are indeed very 
difficult to map due to observational surface brightness 
limitations. 

The radii thi and become smaller towards higher 
redshift. This is a direct consequence of the cosmic evo- 
lution of the virial radii Tvir of the haloes in the Millen- 
nium Simulation, which a ffects the disk sc ale radius fdisk 
in Eq. (fT6)) . As shown bv lMo et"an ()1998D . r^ir scales as 
(1 + z)~^-^ for a fixed circular velocity or as (1 -I- z)~^ 
for a fixed halo mass, consistent with high-redshift ob- 
servations r£^_2;^^^_6]in the Hubble Ultra Deep Field 
(UDF) byl ouwens et al.l ( [20 04) . Their selection criteria 
include all but the reddest starburst galaxies in the UDF 
and some evolved galaxies. It should be emphasized that 
the phenomenological size evolution of galaxies is not 
properly understood, and even modern A^-body/SPH- 
simulation cannot yet accurately reproduce the sizes of 
galaxies. 

In analogy to the mass-radius power-law relation for 
HI, our simulation predicts a similar relation, again 
nearly a square-law, for H2 at redshift z = (see Fig.[7Ji). 
This power-law is consistent w ith observations of th e two 
face-on spi ral galaxies M51 (ISchuster et al.ll2007t l and 
NGC 6946 (jCrosthwaite fc Turner] l2007f) . For the small 
H2/HI-ratios found in the local Universe the A/h2^''H2 
relation is linked to the Mm-rm relation, because both 
Mfjj and rn^ can be regarded as a fraction (< 1) of, 
respectively, Mhi and rui. However, there is no funda- 
mental reason for a constant surface density of H2 and 
the smaller sizes of high-redshift galaxies implies a higher 
pressure of the ISM and thus a much higher molecular 
fraction by virtue of Eqs. (flOl [T5|) . Therefore, H2-masses 
become uncorrelated to HI and tend to increase with 
redshift out to z ~ 5, while decreases. Hence, the 
-Mh2~''H2 relation must move away from the power-law 
found at z = (see Figs. [Hi-f). 

5. REALISTIC VELOCITY PROFILES 

In this section we derive circular velocity profiles and 
atomic and molecular-radio line profiles for the simulated 
galaxies in the DeLucia-catalog. Circular velocity pro- 
files Vc{r) for various galaxies are derived over the Sec- 
tions I5.1H5.3I and transcribed to radio line profiles for 
edge-on galaxies in Section [5T4l Results for the local and 
high-redshift Universe are presented in Section 15.51 



5.1. Velocity profile of a spherical halo 

To account for the narrowness of the emission lines ob- 
served in the central gas regions of many galaxies (e.g. 



iSautv et al.l[2003l) . we require a halo model with van- 
ishing velocity at the center, as opposed to, for exam- 
ple, the commonly adopted singular isothermal sphere 
with a density Phaioi^) ~ f^'^ and a constant veloc- 
ity profile. We chose the Navarro-Frenk- White (NFW, 
iNavarro et al.l Il995l 119961 ) model, which relies on high 
resolution numerical simulations of dark matter haloes 
in equilibrium. These simulations revealed that haloes 
of all masses in a variety of dissipation-less hierarchical 
clustering models are well described by the spherical den- 
sity profile 



(31) 



where po is a normalization factor and rg is the char- 
acteristic scale radius of the halo. This profile is also 
supported by the Hubble Space Telescope analysis of the 
weak lensing induced b y the galaxy cluster MS 2053- 
04 at redshift z = 0.58 (|Hoekstra et al.ll2002[) . /9haio(?-) 
varies as r~^ at the halo center and continuously steep- 
ens to for r 00. It passes through the equi- 
librium profile of the self-gravitating isothermal sphere, 
i.e. Ph2.\o{r) r~2, at r = Tg. 

The definition of rvir in the Millennium Simulation 
given in Eq. (|17p implies that 



Po 
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~3~ MX 
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halo 



(32) 



■ Chalo) — Chalo/ (1 + Chalo ) 

where Chaio = ^vir/fs is referred to as the halo con 
centration parameter. Most numerical models 
diet that Chalo scales with the virial mass Mvir, 



pre- 
de- 



nned as tne mass msi ae trie raams ryjr, accoramg to a 
power- law ( e.g. Navarro et al.lll997t iBullock et al.ll2"00ll : 
iDolag et al] [2004: Hennawi et al."2007^. Here we shaU 
use the result of Hcnnawij;_t_aL (20_07) , 



Chalo 



12.3 



(33) 



1 + z VL3- 10i3/i-iMo, 

which is consistent with recent empirical values of 
the matter concentration in galaxy clusters derived 
from X-ray measurements a nd strong lensing data 
(jComerford fc NataraiMll2007| ). 
For a spherical halo, the circular velocity profile is 

given by V,''-'°\r) = GMhaio(r)/r with Mhaio(r) = 
47r drr^ Phaio{r). Using Eqs. ([?T1[5^ . this implies that 



halo^ 



ix) 



GA'U 



ln(l -I- ChaloX) - 



a;[ln(l -|- Chaio) 



(34) 



l + ChaloJ 

where x = r/rvir (thus Chaio^^ — f/Ta)- This veloc- 
ity vanishes at the halo center, then climbs to a max- 
imal value l^iax = 1.65rsVGpo a.t r — 2.16 Tg, from 
where it decreases monotonically with r, typically reach- 
ing 0.65 — 0.95 Knax at r = Tvir {x = 1) with the extremes 
corresponding, respectively, to Chaio — 25 and Chaio — 5. 
For larger radii, the velocity asymptotically approaches 
the point-mass velocity profile (r) = GM^i^/r. 

5.2. Velocity profile of a flat disk 

For simplicity, we assume in this section that the galac- 
tic disk is described by a single exponential surface den- 
sity for stars and cold gas, 

S^^'W^TT-^exp— , (35) 



fdisk 
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where M'^^^^ is the total disk mass, taken as the sum 
of the cold gas mass and the stellar mass in the disk. 
In most real galaxies the stellar surface densities, are 
slightly more compact (see Section 13.11) , but we found 
that including this effect docs not significantly modify 
the shape of the atomic and molecular emission lines. 
In fact, the radius, which maximally contributes to the 
disk mass, i.e. the maximum of r E'^"''^(r), is r = rdisk- 
Therefore, we expect the gravitational potential to differ 
significantly from the point-mass potential only for r of 
order rdisk or smaller. Applying Poisson's equation to the 
surface density of Eq. (|35|) . the gravitational potential in 
the plane of the disk becomes 

^ ^_ GM'^'^^ ff exp(-f/rdisk)rdrd(^ 
'^disklrj- 2^r2.^, (r2+r~2-2rr~ cos 0)1/2' ^"^^^ 

where the integration surface D is given by f G [0, cxd), 
e = [0,27r). 

The velocity profile for circular orbits in the plane of 

the disk can be calculated as V^,'^'^'^ = r d(pdisk/dr. The 
integral in Eq. (j36[) is elliptic, and hence there are no 
exact closed- from expressions for tpdisk and I^'^'"'^. How- 
ever, in this study we numerically found that an excellent 
approximation is given by 

g^^disk 



V 



disk^ 



{x) 



Cdisk + 4.8cdiskexp[-0.35ediska:^ - 3.5/(cdiska:)] 
Cdisk2^ + (cdiska;)"^ + 2(cdiska;)~^/^ 



(37) 
(38) 



where Cdisk = ''vir/^disk is the disk concentration param- 
eter in analogy to the halo concentration parameter Chaio 
of Section [Ql Eq. ([37|) is accurate to less than 1% over 
the whole range r = 0— 10 rdisk and it correctly converges 
towards the circular velocity of a point-mass potential, 

l/h'^i°^(r) = GM^'^'^/r, for r ^ oo. 

Like in Section 13.21 we shall use Eq. (fT6|) to compute 
rdisk and Cdisk in the DeLucia-catalog. This approach 
is slightly incon sistent because Eq. (flB)) was derived by 
iMo et ail (|1998 f) under the assumption that the disk is 
supported by an isothermal halo with /9haio('') ^ f '^, 
while in Section [5.11 we have assumed the more complex 
NFW-profile. We argue, however, that Eq. (fT6|) with the 
empirical correction of Eq. ((20|) is sufficiently accurate, as 
it successfully reproduces the observed relation between 
-^bary and rdisk (see Fig. [2]) as well as the relation be- 
tween Afni and rni (see Fig. [H). It can also be shown 
that, for realistic values of the halo concentration Chaio 
(10 — 20 for one- galaxy systems) and the spin parameter 
A (0.05-0.1, e.g. lMo et anil998f) . the scale radius of the 
halo rs and the disk radius rdisk are similar. Hence, the 
main mass contribution of the disk comes indeed from 
galactocentric radii, where the halo profile is approxi- 
m ately isoth e rmal, thus justifying the assumption made 
bv lMo et al.] (jl998[ ) to derive Eq. (fTe]) . 

5.3. Velocity profile of the bulge 

Many models for the surface brightness or surface den- 
sit y profiles of bulges h ave been proposed (e.g. overview 
bv lBalcells erall l200l. A rough consensus seems es- 
tablished that no single surface density profile can de- 
scribe a majority of the observed bulges, but that 
they are generally well matched by the class of Sersic- 
functions llSersidll968fl . Eb"ig<=(r) - exp[-(r/rbuigc)^/"], 



where the exponent n depen ds on the morphological 
type (jAndredakis et al.l Il995f ). such that n « 4 for 
lenticular/early- type galaxies (de Vaucouleurs-profile) 
and n « 1 fo r the bulges of late-typ e galaxies (exponen- 
tial profile). ICourteau et al.l (1996) find slightly steeper 
profiles with n = 1 — 2 for nearly all spirals in a sample of 
326 spiral galaxies using deep optical and IR photometry, 
and they show that by imposing n = 1 for all late-type 
galaxies, the ratio between the exponential scale radius 
of the bulge rbuigc and the scale radius of the disk rdisk 
is roughly constant, rbuigo ~ O.lfdisk ~ 0.05 rdisk- We 
shall therefore assume that all bulges have an exponen- 
tial projected surface density. 



^bulgo 
2^^bulgc 



exp 



rbulge 



(39) 



with rbuigc = 0.05 rdisk- 

For simplicity, we assume that the bulges of all galaxies 
are spherical and thus described by a radial space den- 
sity function phuige{f). This function is linked to the pro- 
jected surface density via S^'^'sc^,^) = dz pbuigo [(r^-l- 

z^)i/^]. Numerically, we find that this model for 
Pbuige(^) is close ly approximated by the Plummcr model 
(|Plunimeilll911f) . more often used in the context of clus- 
ters, 



Pbulgo(r) 



3 M^"'so 

^'^ '^Plummcr 



rpi 



ummcr 



-1 -5/2 



(40) 



lummer 



with a characteristic Plummer radius rpiu 
1.7 rbuigc. The circular velocity profile V^^^^° correspond- 
ing to Eq. (HOI) is given by T/buig<=^(r) = GM^''^<i%r)/r 
with A/'^^'s'=(r) = An df pbuigo(r). This solves to 



bulgo2/ N _ GA/'^^'g" (Cbulgea:)^Cbulgc 
'^vir [1 + {C^ulso^rf 



(41) 



where Cbulgo = ?'vii /?'Plummcr 

centration parameter 



12 Cdisk is the bulge con- 



5.4. Line shapes from circular velocities 

The addition rule for gravitational potentials implies 
that the circular velocity profile of the combined halo- 
disk-bulge system in the plane of the disk is given by 



V^'ix) = ix) + Vr'' (x) + T/r'sc (^)^ (42) 

where x = r/rvir as in Sections I5.1H5.3I According to 
Eqs. (|341 [571 HT]) . this profile is determined by six param- 
eters: the three form-parameters Chaio, Cdisk, Cbuigo and 
the three mass-scales M^-„/r^-„, M'^'''^ / r^-„ , M'^^'s^/rvir. 
The form-parameters were calculated as explained in Sec- 
tions l5.lH5.3l while the mass-scales were directly adopted 
from the DeLucia-catalog. For the satellite galaxies with 
no resolved halo (see Section [2]), A/vir and rvir were ap- 
proximated as the corresponding quantities of the origi- 
nal galaxy halo just before its disappearance. An exem- 
plar circular velocity profile for a galaxy in the DeLucia- 
catalog at redshift z = is shown in Fig. [51 

In order to evaluate the profile of a radio emission 
line associated with any velocity profile Vc{r), we shall 
first consider the line profile of a homogeneous flat ring 
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Fig. 9. — Circular velocity profile of a typical simulated galaxy 
with a small bulge at redshift 2 = 0. The total circular velocity 
(solid line) is given by the circular velocity of the halo (dashed 
line) , th e disk (dash-dotted line), and the bulge (dotted line) via 
Eq. gg. 




Fig. 10. — Apparent velocity Vobs induced by the infinitesimal 
ring element at the angle 7. 



with constant circular velocity Vc and a total luminos- 
ity of unity. If a point of the ring is labeled by the 
angle 7 it forms with the line-of-sight (see Fig. [rU]) . 
the apparent projected velocity of that point is given 
by Kbs = K sin 7. The ensemble of all angles 7 € 
[0, 27r) therefore spans a continuum of apparent veloci- 
ties I4bs G (— with a luminosity density distribu- 
tion -iACKibs, K;) ~ d7/dVi3bs- Imposing the normalization 
condition / dV^bs'0(K)bs) 1, we find that the edge-on 
line profile of the ring is given by 







if iKbsl < V. 
otherwise. 



(43) 



This profile exhibits spurious divergent singularities at 
|K)bs| which, in reality, are smoothed by the ran- 

dom, e.g. turbulent, motion of the gas. We assume that 
this velocity dispersion is given by the constant Cgas = 
8 km s~^, which is consistent with the velocity disper- 
sions observBd_a£ross_yie_disks of several nearby galaxies 
(e.g. IShostak fc van der Kruitl[l98l iDickev etatl 119901 : 
iBurtonI Il971h . The smoothed velocity profile is then 



V'(Kbs,K) = 



dV^exp 



(Kbs-^)' 
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^{V,V,), (44) 



which conserves the normalization / dKibsV'l^bs) = 
1. Some examples of the functions ■(/'(^Dbs, K) and 
■0(V"obs, Vc) are plotted in Fig. [TT] 
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Fig. 11. — Illustration of the functions ip (Eq.\^ and i/i (Eq.liSt. 
which represent the normalized emission line of a homogeneous 
edge-on disk or ring with constant circular velocity. 

From the edge-on line profile V'lKabsj K) of a single ring 
and the face-on surface densities of atomic and molecular 
gas, SHi(r) and (?"), we can now evaluate the edge-on 
profiles of emission lines associated with the entire HI- 
and Il2-disks, respectively. Since H2-densities are most 
commonly inferred from CO-detections, we shall here- 
after refer to all molecular emission lines as "the CO- 
line" . The edge-on line profiles (or "normalized luminos- 
ity densities") *Hi(K>bs) and *co(Kbs) are given by 



*Hl(K,bs) = 



277 



drrSHi(r)V(K,bs,K(r-)), (45) 

JWhi Jo 
27r 

*Co(K,bs) = T;j— / drrSH,(r)V(Kb.,T4(r)).(46) 

These two functions satisfy the normalization conditions 
/ dyobs*Hi(14bs) = 1 and / dyob.«'co(K>b«) = 1- To 
obtain intrinsic luminosity densities, ^'Hi(KDbs) must be 
multiplied by the integrated luminosity of the Hl-line [see 
Eq. (f2T|) ] and ^'co(K)bs) must be multiplied by the in- 
tegrated luminosity of the considered molecular emission 
line, e.g. the integrated luminosity of the CO(l-0)-line 
given in Eq. ((23|) . 

Fig. [121 displays the line profiles 5'Hi(K)bs) and 
^'co(K)bs) for the exemplar galaxy with the velocity 
profile shown in Fig. [9] All line profiles produced by 
our model are mirror-symmetric, but they can, in prin- 
ciple, differ significantly from the basic double-horned 
function '4'{Vohs)- Especially for CO, where the emission 
from the bulge can play an important role, several lo- 
cal maxima can sometimes be found in the line profile, 
in qualitative agreemen t with various observations (e.g. 
iLavezzi fc Dickevlll998l ). 

5.5. Results and discussion 

For every galaxy in the DeLucia-catalog, we computed 
the edge-on line profiles ^'Hi(Kbs) and 4*00 (Kbs), from 
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Fig. 12. — Simulated edge-on HI- and CO-emission lines for the 
exemplar galaxy, for which the circular velocity profile i s sh own in 
Fig. [9] The line profiles have been computed using Eq. (1451 . 

which we extracted the line parameters indicated in 
Fig.ini ^-^j EE *Hi(0) and = *co(0) are the lumi- 
nosity densities at the hne center, and and ^'co'^ 
are the peak luminosity densities, i.e. the absolute max- 
ima of ^'Hi(Vobs) and *co(Kbs)- w^i""^ and w^°q^ are 
the hne widths measured between the left and the right 
maximum. These values vanish if the line maxima are 
at the line center, such as found, for example, in slowly 
rotating systems, w^^, Wqq, w^I, and w^q are the line 
widths measured at, respectively, the 50-percentile level 
or the 20-percentile level of the peak luminosity densi- 
ties - the two most common definitions of observed line 
widths. 

We shall now check the simulated line widths against 
observations by analyzing their relation to the mass of 
the galaxies. Here, we shall refer to all line width versus 
mass relations as Tully-Fisher relations (TFRs), since 
they are generalized versions of the original relation be- 
twe en line widths and opt ical magnitudes of spiral galax- 
ies ()Tullv fc Fished I1977D . A variety of empirical TFRs 
have been publishe d, such as the ste llar mass-TFR and 
the baryonic-TFR (jMcGaugh et al.l i2000). The latter 
relates the baryon mass (stars-|-gas) of spiral disks to 
their line widths (or circular velocities) and is probably 
the most fundamental TFR detected so far, obeying a 
single power-law over five orders of magnitude in mass. 
We have also investigated the less fundamental empirical 
TFR between Mm and - hereafter the HI-TFR - us- 
ing the spiral galaxies of the HIPASS catalog. Assuming 
no prior knowledge on the inclinations of the HIPASS- 
galaxies, but taking an isotropic distribution of their axes 
as given, we found that the most-likely relation is 

for the Hubble parameter h = 0.73. Relative to Eq. (|47p 
the HIPASS data exhibit a Gaussian scatter with a = 
0.38 in log(MHi). Our method to find Eq. ^ will be 
detailed in a forthcoming paper (Obreschkow et al. in 
prep.), especially dedicated to the Hl-TFR. 
Figs.fT3b-d show four TFRs at redshift z = 0. Each fig- 



ure represents 10"^ simulated galaxies (black), randomly 
drawn from the simulation with a probability propor- 
tional to their cold gas mass in order to include the rare 
objects in the high end of the MF. Spiral and elliptical 
galaxies are distinguished as black dots and open circles. 

Fig.fTSk shows the simulated Hl-TFR together with the 
empirical counterpart given in Eq. (|17|) . This compari- 
son reveals good consistency between observation and 
simulation for spiral galaxies. However, the elliptical 
galaxies lie far off the Hl-TFR. In fact, simulated el- 
liptical galaxies generally have a significant fraction of 
their cold hydrogen in the molecular phase, consistent 
with the gala xy-type dependen ce of the H2 /Hl-ratio first 
identified bv lYoung fc Knezekl |l989). Therefore, HI is 
a poor mass tracer for elliptical galaxies, both in simu- 
lations and observations, leading to their offset from the 
TFR when only Hl-masses are considered. There seems 
to be no direct analog to the HI-TFR for elliptical galax- 
ies. 

Figs. [T5b . c respectively display the simulated stel- 
lar mass-TFR an d the baryonic-TFR, tog ether with the 
observed data of iMcGaugh et al.l ([2000) corrected for 
h — 0.73. These data include various galaxies from 
dwarfs to giant spirals, whose edge-on line widths were 
estimated from the observed ones using the inclinations 
determined from the optical axis ratios. Figs. 113b . c 
reveal a surprising consistency between simulation and 
observation. In Fig. [T5b . both the simulated and ob- 
served data show a systematic offset from the power-law 
relation for all galaxies with < 200 kms~^. Yet, 
the power-law relation is restored as soon as the cold gas 
mass is added to the stellar mass (Fig. [T3b). thus con- 
firming that the TFR is indeed fundamentally a relation 
between mass and circular velocity. 

It is interesting to consider the prediction of the simula- 
tion for the most fundamental TFR, i.e. the one between 
the total dynamical mass, taken as the virial mass A/vir i 
and the circular velocity, represented by the line width 
This relation is shown in Fig. [T3H and reveals in- 
deed a 2-3 times smaller scatter in log(mass) than the 
baryonic-TFR, hence confirming its fundamental charac- 
ter. 

Although the simulated elliptical galaxies shown in 
Figs. [T5b-d roughly align with the respective TFRs for 
spiral galaxies, their scatter is larger. This is caused by 
the mass-domination of the bulge, which leads to steep 
circular velocity profiles Vc{r) with a poorly defined ter- 
minal velocity. Therefore, line widths obtained by aver- 
aging over the whole elliptical galaxy are weak tracers of 
its spin. This picture seems to correspond to observed 
elliptical galaxies, where the central line widths, corre- 
sponding to the velocity dispersion in the bulge domi- 
nated parts, are more correlated to the stellar mass than 
the line wid ths of the whole galaxy (see Faber- Jackson 
relati on, e.g. iFaber fc Jacksonlll976l ) . 

KassineLaU (12Q07D no"ted that S'0.5 = {0.5V^ + 
cLo)^^^ is a better kinematic estimator than the circu- 

gas / 

lar velocity Vc, in the sense that it markedly reduces 
the scatter in the stellar mass-TFR. However, since our 
model assumes a constant gas velocity dispersion ag^s for 
all galaxies, it is not possible to investigate this estima- 
tor. 

The predicted evolution of the four TFRs in Figs.fTSk- 
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Fig. 13. — Relations between edge-on line widths and different mass tracers for the local Universe. 10'^ simulated galaxies are represented 
by black dots (spiral galaxies) and black circles (elliptical galaxies). The solid lines represent power-law fits to the simulated spiral galaxies; 
their slopes are respectively qhi = 2.5, Ostars = 4.3, abary = 3.9, and avir = 3.6. In case of Fig. 113b . this fit only includes galaxies with 
Mstars > 10^ Mq. Fig. I13H does not include satellite galaxies without haloes (see Section [2J, for which Afvir is poorly defined. The dashed 
red line and shaded zone in Fig. 113b , represent our obs erva tional determination and l-cr scatter of the HI-TFR from the HIPAS S data (see 
Sectio n |5.5| I. The rose dots and dashed lines in Figs. [TSb . c are the observational data and power-law regressions from McGaugh et al.l 
((2003) and referen ces t herein; this sample include low surface brightness galaxies. The slopes of these empirical power-laws are respectively 
OHI = 2.8 (see Eq.|47J, Ostars = 4.4, and abary = 4.4. 



d is shown in Figs.[T4k-d. In all four cases, the simulation 
predicts two important features: (i) galaxies of identi- 
cal mass (respectively Mhi, Mgtars, Afbary, -^^vir) have 
broader lines (and larger circular velocities) at higher 
redshift, and (ii) the scatter of the TFRs generally in- 
creases with redshift. The first feature is mainly a 
consequence of the mass-radius-velocity relation of the 
dark matte r haloes assumed i n the Millennium Simu- 
lation (see ICroton et al] l2006t iMo et al.l |1998[ ). This 
relation predicts that, given a constant halo mass, Vc 
scales as (1 -I- zy^^ for large z. Furthermore, the ratios 
Mm/Mvir and Mstars/Mvir on average decrease with in- 
creasing redshift, explaining the stronger evolution found 
in Figs. [T4k. b relative to Figs. [T4b. d. 

The increase of scatter in the TFRs with redshift is a 
consequence of the lower degree of virialization at higher 
redshifts, which, in the model, is accounted for via the 
spin parameter A of the haloes. A is more scattered at 
high redshift due to the young age of the haloes and the 
higher merger rates. More scatter in A leads to more 
scatter in the radius rdisk via Eq. (jl6p and thus to more 
scatter in the circular velocity Vc via Eqs. ([37l |4T |) . 

Current observational databases of resolved CO-line 
profiles are much smaller than Hl-databases and their 
signal/noise characteristics are inferior. Nevertheless ef- 
forts to_ch£ck_tiwuserf widths for probing TFRs 
(e.g. iLavezzi fc Dickey! 1 1998( ) have led to the conclusion 
that in most spiral galaxies the CO-line widths are very 
similar to Hl-line widths, even though the actual line 
profiles may radically differ. Fig. [15] shows our simu- 
lated relation between and Wqq, as well as the lin- 
ear fit to observations o f 44 galaxies in different clusters 
(jLavezzi fc Dickeylll998( ). These observations are indeed 
consistent with the simulation. The simulated ellipti- 
cal galaxies tend to have higher w^Q/w^^-ratios than 
the spiral ones, due to fast circular velocity of the bulge 
component implied by its own mass. 



The line profiles and widths studied in this section cor- 
respond to galaxies observed edge-on. First order correc- 
tions for spiral galaxies seen at an inclination i ^ 90 deg 
can be obtained by dividing the normalized luminosity 
densities ^^j, "^^f^, ^co^ ^co" ^y sini, while multiply- 
ing the fine widths w^J}, w'^^, w^r"^, w^q, w'^q, w^q^ 
by sinj. More elaborate corrections, accounting for the 
isotropy of the veloc ity dispersi on CTgas, we r e give n by 
iMakarov etall ()1997D and Kann appan et al.l (|2002[ ). 

6. DISCUSSION 

We used a list of physical prescriptions to post-process 
the DeLucia-catalog and showed that many simulation 
results match the empirical findings from the local Uni- 
verse. However, this approach raises two major ques- 
tions: (i) Are the applied prescriptions consistent with 
the DeLucia-catalog in the sense that they represent a 
co mpatible extension o f the semi- analytic recip es used 
by lDe Lucia fc BlaizotI (|2007D and lCroton et all (|2006D ? 
(ii) What are the limitations of our prescriptions at low 
and high redshifts? 

6.1. Consistency of the model 

The consistency question arises, because the DeLucia- 
catalog relies on a simplified version of a Schmidt- 
Kennicutt law (|Schmidtll 19591 : lKennicuttl [l998l). I.e. a pre- 
scription where the star formation rate (SFR) scales as 
some power of the surface density of the ISM. However, in 
a smaller-scaled picture, new stars are bred inside molec- 
ular clouds, and hence it must be verified whether our 
prescription to assign H2-masses to galaxies is compat- 
ible with the macroscopic Schmidt-Kennicutt law. Our 
prescription exploited the empirical power-law between 
the pressure of the ISM and its mol ecular content, as first 
presented by ,Bhtz fc Rosolowskvl (l2b04. ,2006). Based 
on this power-law relation, iBlitz fc Rosolowskvl (|2006f ) 
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Fig. 14. — Simulated cosmic evolution of the different line width-mass relations shown in Fig. 1131 Sp iral and elliptical galaxies are 
respectively represented by dots and circles. Black color corresponds to redshift z = (identically to Fig. I13|l . while blue and red color 
respectively represent z = 4.89 and z = 10.07. The solid black lines are power-law fits to the spiral galaxies at ^ = 0, where in case of 
Fig. 114b only galaxies with Mstars > 10^ Mq have been considered. The number of elliptical galaxies decreases with redshift — a consequence 
of the merger- and instability-driven prescriptions for bulge formation in the DeLucia-catalog. 
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Fig. 15. — Relation between line widths of HI and CO. lO'^ 
simulated galaxies are represented by black dots (spiral galax- 
ies) and black circles (elliptical galaxies). The red dashed line 
and rose-shaded zone represent the best fit and its l-cr confidence 
interval to o bservations of 44 galaxies in clusters presented by 
ILavezzi fc Di ckey (1998). 

themselves formulated an alternative model for the com- 
putation of SFRs in galaxies, which seems more funda- 
mental than the Schmidt-Kennicutt law. Applying both 
models for star formation to six molecule-rich galaxies in 
the local Universe, they showed that their new pressure- 
based law predicts SFRs nearly identical to the ones pre- 
dicted by the Schmidt-Kennicutt law. Therefore, our 
choice to divide cold hydrogen in HI and H2 accord- 
ing to pressure is inde ed consistent with t h e pre scrip- 
tion for SFRs used by iDe Lucia fc Blaizoti (|2007[) and 
ICroton et all (|2006h . 

6.2. Accuracy and limitations at z ^ 



A first limitation of our simulation comes from the as- 
sumption that the surface densities of HI and H2 are 
axially symmetric (no spiral structures, no central bars, 
no warps, no satellite structures). In general, our model 
describes all galaxies as regular galaxies - as do all semi- 
analytic models for the Millennium Simulation. Hence, 
the simulation results cannot be used to predict the HI- 
and H2-properties of irregular galaxies. 

While our models allowed us to reproduce the observed 
relation between Mhi and rni remarkably well for various 
spiral galaxies (e.g. Fig. [5]), it tends to underestimate the 
size of Hl-distributi ons in elliptical galaxie s. For exam- 
ple observations by iMorganti et al.l (|1997f ) show that 7 
nearby E- and SO-type galaxies all have very complex 
Hl-distributions, often reaching far beyond the corre- 
sponding radius of a mass-equivalent disk galaxies. The 
patchy Hl-distributions found around elliptical galaxies 
are probably due to mergers and tidal interactions, which 
could not be modeled in any of the semi-analytic schemes 
for the Millennium Simulation. 

Another limitation arises from neglecting stellar bulges 
as an add itional source of disk-pressure in Eq. (O 
(|Elmeg rccn 1989). Especially the heavier bulges of early- 
type spiral galaxies could introduce a positive correction 
of the central pressure and hence increase the molecu- 
lar fraction, thus leading to very sharp H2-peaks in the 
galaxy centers, such as observed, f or example, in th e 
SBb-type spiral galaxy NGC 3351 (iLerov et al.l 1200^. 
Our model for the H2-surface density of Eq. (fT2)) fails 
at predicting such sharp peaks, although the predicted 
total HI- and H2-mass and the corresponding radii and 
line profiles are not significantly affected by this effect. 

6.3. Accuracy and limitations at z > 

Additional limitations are likely to occur at higher red- 
shifts, where our models make a number of assumptions 
based on low-redshift observations. Furthermore, the un- 
derlying DeLucia-catalog itself may suffer from inaccura- 
cies at high redshift, but we shall restrict this discussion 
to possible issues associated with the models in this pa- 
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per. 

Regarding the subdivision of hydrogen into atomic 
and molecular material (Section [3]), our most critical as- 
sumption is the treatment of all galactic disks as regu- 
lar exponential structures in hydro-gravitational equilib- 
rium. This model is very likely to deviate more from 
the reality at high redshift, where galaxies were gener- 
ally less virialized and mergers were much more frequent 
(|de Ravel et al.l [20081 ). Less virialized disks are thicker, 
which would decrease the average pressure and fraction of 
molecules compared to our model. Yet, galaxy mergers 
counteract this effect by creating complex shapes with 
locally increased pressures, where H2 forms more effi- 
ciently, giving rise to merger-driven starbursts. There- 
fore, it is unclear whether the assumption of regular disks 
tends to underestimate or overestimate the H2/HI-ratios. 

Another critical assumption is the high-rcdshift valid- 
ity of the local relation between the H2/HI-ratio i?moi 
and the external gas pressure P (Eq. [6|) . This relation 
is not a fundamental thermodynamic relation, but rep- 
resents the effective relation between the average H2 /HI- 
ratio and P, resulting from complex physical processes 
like cloud formation, H2-formation on metallic grains, 
and H2-destruction by the photodissociative radiation 
field of stars and supernovae. Therefore, the i?moi — P 
relation could be subjected to a cosmic evolution result- 
ing from the cosmic evolution of the cold gas metallicity 
or the cosmic evolution of the photodissociative radia- 
tion field. However, the metallicity evolution is likely to 
be problem only at the highest redshifts (z > 10). Ob- 
servations in the local Universe show that spiral galax- 
ies with metallicities diffe ring by a factor 5 fall o n the 
same i?moi — P relation (jBlitz fc Rosolows"kvll2006f ) . Yet, 
the average cold gas metallicity of the galaxies in the 
DeLucia-catalog is only a factor 1.9 (3.6) smaller at 2; = 5 
(z = 10) than in the local Universe. These predictions 
are consistent with observational evidence from the Sloan 
Digital Sky Survey (SDSS) that stellar metallicities were 
at most a factor 1.5-2 smaller at z 3 than today 
(jPanter et al.|[200l . The effect of the cosmic evolution of 
the photodissociative radiati on field on the i?moi — P re - 
lation is difficult to assess. iBlitz fc Rosolowskvl ()2006f ) 
argued that the ISM pressure and the radiation field 
both correlate with the surface density of stars and gas, 
and therefore the radiation field is correlated to pressure. 
This is supported by observations in the local Universe 
showing that the i?,noi — P relation holds true for dwarf 
galaxies and spiral galaxies spanning almost three orders 
of magnitude in SFR. For those reasons, the i?moi — P 
relation could indeed extend surprisingly well to high red- 
shifts. 

In the expressi on for the disk-pressure in Eq. ([7]) 
(jElmegreenl ll989D . we assumed a constant average 
velocity dispersion ratio {fcr)- Observations sug- 
gest that Vr/a^Fis decre ases si gnificantly w ith redshift 
(jForster Schreiber et all 120061: ICenzel et all [2008^. and 
therefore (/o-) perhaps increases. This would lead to even 
higher H2/HI-ratios than predicted by our model. How- 
ever, according to Eq. (fTU]) this is likely to be a problem 
only for galaxies with Mgtars > -^gas, while most galaxies 
in the simulation at z > 2 are indeed gas dominated. 

Regarding cold gas geometries and velocity profiles, the 
most important limitation of our model again arises from 
the simplistic treatment of galactic disks as virialized ex- 



ponential structures. Very young galaxies (< 10^ yrs) or 
galaxies undergoing a merger do not conform with this 
model, and therefore the predicted velocity profiles may 
be unreal and the disk radii may be meaningless. This 
is not just a limitation of the simulation, but it reveals a 
principal difficulty to describe galaxy populations domi- 
nated by very young or merging objects with quantities 
such as rni or which are common and useful for 

isolated systems in the local Universe. 

The radio line widths predicted by our model (Sec- 
tioning!]) may be underestimated at high redshift, due to 
the as sumption of a constant r a ndom velocity dispersion 
cTgas- iForster Schreiber et al.l (|2006l ) found 14/CTgas ~ 
2 — 4 for 14 UV-selected galaxies at redshift z w 2. 
This result suggests that radio lines at z « 2 should 
be about 20%-30% broader than predicted by our model, 
and therefore the evolution of the TFRs could be slightly 
stronger than shown in Fig. [141 

In summary, the HI- and H2-properties predicted for 
galaxies at high redshift are generally uncertain, even 
though no significant, systematic trend of the model- 
errors could be identified. Perhaps the largest deviations 
from the real Universe occur for very young galaxies or 
merging objects, while isolated field galaxies, typically 
late-type spiral systems, might be well described by the 
model at all redshifts. 

In Section 13.41 we ascribed CO(l-0)-luminosities to 
the H2-masses using the metallicity dependent X-factor 
of Eq. (|26p . This model neglects several important as- 
pects: (i) the projected overlap of molecular clouds, 
which is negligible in the local Universe, may become sig- 
nificant at high redshifts, where galaxies are denser and 
richer in molecules; (ii) the temperature of the cosmic 
microwave background (CMB) increases with redshift, 
hence changing the level population of the CO-molecule 
(Silk & Spaans 1997; Combes et al. 1999); (iii) the CMB 
presents a background against which CO-sources are de- 
tect; (iv) the molecular material in the very dense galax- 
ies, such as Ultra Luminous Infrared Galaxies, may be 
distributed smoothly rather than in clouds and clumps 
(jPowncs ct al. 1993); (v) the higher SFRs in early galax- 
ies probably led to higher gas t emperatures, hence ch ang- 
ing the CO-level population^. ICombes et al.l ()1999f ) pre- 
sented a simplistic simulation of the cosmic evolution of 
X, taking the cosmic evolution of metallicity and points 
(i) and (ii) into account. They found that for an H2-rich 
disk galaxy {X) increases by a factor 1.8 from redshift 
z = to z = 5. This value approximately matches the 
average increase of X by a factor 2 predicted by our 
simulation using the purely metallicity-based model of 
Eq. (|26l) . This indicates that the effects of (i) and (ii) 
approximately balance each other. If a more elaborate 
model for X becomes available, the latter can be directly 
applied to correct our CO-predictions. In fact, the X- 
factor only affects the CO(l-0)-luminosity Leo calcu- 
lated via Eq. ([23)) . but has no effect on the line properties 

considered in this paper, namely the line widths Wq'q^, 
Wqq, Wqq and the normalized luminosity densities '^qq, 

iTfinax 

This list is not exhaustive, see IMalonev fc Bladd iITqSSI' ) : IWalil 
1I2OO7I) for an overview of the physical complexity behind the X- 
factor. 
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7. CONCLUSION 

In this paper, we presented the first attempt to in- 
corporate detailed cold gas properties in a semi-analytic 
simulation of galaxies in a large cosmological volume. To 
this end, we introduced a series of physical prescriptions 
to evaluate relevant properties of HI and H2 in simulated 
model-galaxies. 

When applied to the DeLucia-galaxy catalog for the 
Millennium Simulation, our recipes introduce only one 
free parameter in addition to the 9 free parameters of 
the semi -analytic model of t he DeLucia-catalog (see Ta- 
ble 1 in ICroton et al.ll2006[ ). This additional parame- 
ter, i.e. the cold gas correction factor C (Section [2]), was 
tuned to the cosmic space density of cold gas in the lo- 
cal Universe. The additional parameter ^, describing the 
transfer of angular momentum from the halo to the disk 
(Section 13. 2p . is not a free parameter for the hydrogen 
simulation, since it is fixed by the baryon mass-scale 
radius relation. In fact, we deliberately did not adjust 
£ to match the observed HI-mass-HI-radius relation of 
IVerheiienI (|2001l ). in order to check the reliability of our 
models against this observation. 

Based on the DeLucia-catalog, we produced a virtual 
catalog of ~ 3 • 10^ per redshift-snapshot with various 
cold gas properties. This catalog represents an extension 
of the DeLucia-catalog, and it can be used to investi- 
gate a broad variety of questions related to HI, H2, CO 
and their cosmic evolution. The results presented in this 
paper have been restricted to some important examples, 
most of which could be compared directly to available 
observations and hence constitute key results for the ver- 
ification of our simulation: 

1. Based on a pressure-model for the molecular con- 
tent of cold gas, the simulation simultaneously re- 
produces the HI-MF and the H2-MF (resp. the CO- 
LF) observed in the local Universe within the mea- 
surement uncertainties (Fig. [3|) . 

2. The simulated HI-MFs for spiral and elliptical 
galaxies considered individually also match the ob- 
servations for simulated galaxies with well-defined 
galaxy types (Fig.[5|). 

3. The simulated Hl-radii, imply a mass-radius rela- 
tion for HI that matches the empirical counterpart 
(Fig. [8]) , thus confirming that the relation between 
A/hi and rpji is such that the average Hl-density 
inside rni is 3.8 Mq pc~^ for all galaxies in the lo- 
cal Universe, although this value sensibly depends 
on the definition of rm • 

4. The simulation predicts that the mass-radius re- 
lations for HI and H2 are similar in the local Uni- 
verse, but that their high-redshift evolution is com- 
pletely different (Fig. [7]). 



5. The simulated widths of the Hl-radio emission lines 
of spiral galaxies are consistent with the empirical 
HLTFR derived from the HIPASS spiral galaxies 
(Fig. [T3k): and the simulation predicts that there 
is no analog HI-TFR for elliptical galaxies. 

6. The simulated the stellar mass-TFR and the bary- 
onic TFR reveal good agreement with the empirical 
TFRs for both spiral and elliptical galaxies in the 
simulation fFigs. [T5b. c). 

7. These TFRs are observable manifestations of a 
more fundamental relation between circular veloc- 
ity and total dynamical mass, as suggested by the 
small scatter in the relation between and Mvir 
(Fig. [I3H). 

8. At higher redshift, the simulation predicts that the 
above TFRs remain valid (except for HI at z w 10), 
but that their scatter increases and their zero-point 
is shifted towards higher velocities at fixed mass 
- a fundamental prediction of hierarchical growth 
(Fig. [HI). 

The good match between simulation and observation 
regarding gas masses, disk sizes, and velocity profiles 
supports the models and recipes established in this pa- 
per. It also validat e s the semi-analytic recipes used by 
iDe Lucia fc BlaizoH (|2007D and s upports the Millennium 
Simulation ("Sp ringel et al.l 120051 ) as a whole. 

In forthcoming investigations the presented extension 
of the DeLucia-catalog towards cold gas properties could 
be used to investigate more elaborate questions. For 
example, what is the bias of the cosmic structure, for 
example of the power spectrum, revealed in Hl-surveys 
or CO-surveys compared to the underlying dark mat- 
ter structure? How many Hl-sources can we expect to 
detect in future experiments performed by the SKA? Or 
how does the global H2/HI-ratio evolve with redshift and 
how does it relate to the observed evolution of the SFR 
density? 
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